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Abstract. We study mapping properties of operators with kernels defined via 
a combination of continuous and discrete orthogonal polynomials, which pro- 
vide an abstract formulation of quantum (<?-) Fourier type systems. We prove 
Ismail's conjecture regarding the existence of a reproducing kernel structure 
behind these kernels, by establishing a link with Saitoh 's theory of linear trans- 
formations in Hilbert space. The results are illustrated with Fourier kernels 
with ultraspherical weights, their continuous g-extensions and generalizations. 
As a byproduct of this approach, a new class of sampling theorems is obtained, 
as well as Neumann type expansions in Bessel and 5-Bessel functions. 



1. Introduction 

The Gegenbauer expansion of the two variable complex exponential in terms 
of the ultraspherical polynomials 



(1.1) e" 4 = I» ( I j ]T i\v + k)J u+k {t)Ct{x) 

has the remarkable feature of being at the same time an expansion in a Neumann 
series of Bessel functions. The usefulness of this expansion was made very clear in 
a paper authored by Ismail and Zhang, where it was used to solve the eigenvalue 
problem for the left inverse of the differential operator, on L 2 spaces with ultras- 
pherical weights [20] . The consideration of the g-analogue of this diagonalization 
problem led the authors to extend Gegenbauer 's formula to the g-case. This task 
required the introduction of a new q-analogue of the exponential, a two variable 
function denoted by £ q (x;t) which became known in the literature as the curly q- 
exponential function, bearing the name from its notational convention. Ismail and 
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Zhang's formula is 
(1.2) 

The functions involved in this formula will be defined in section 4. Since its intro- 
duction, the function £ q was welcomed as a proper q-analogue of the exponential 
function, since it was suitable to provide a satisfactory q-analogue of the Fourier 
theory of integral transformations and series developments. This suitability was 
made concrete by Bustoz and Suslov in [5], where the authors introduced the sub- 
ject of g-Fourier series. Some of the subsequent research activity has been already 
collected in a book [29] . Among recent developments not yet included in this 
book, we quote the orthogonality relations for sums of curly exponential functions 
[23], obtained using spectral methods, and the construction of a g-analogue of the 
Whittaker-Shannon-Kotel'nikov sampling theorem [19] . The designation "Quan- 
tum" has appear often in recent literature on g-analysis, as in the monographs 
[22j and [21] . This designation is very convenient, since g-special functions are 
intimately connected with representations of quantum groups [6] . 

An abstract formulation designed to capture the essential properties of q- 
Fourier type systems was proposed in |15j and we proceed to describe it here. Let 
ll be a measure on the real line and {p n (x)} a complete orthonormal system in 
L 2 (li). Let (tj) be a sequence of points on the real line and assume that {r n (x)} is 
a discrete orthonormal system with orthogonality relation 



3=0 

and with dual orthogonality 

oo 

^Vfc(i„)rfc(i m ) = 



fe=0 



P(tn) 



Assume also that the the system {r n (x)} is complete in l 2 (J2 p(tj)Stj)- Now define 
a sequence of functions {F n (x)} by 

oo 

(1.3) F n (x) =J2^k{t n )p k 

k=0 

where {u k } is an arbitrary sequence of complex numbers in the unit circle. The fol- 
lowing theorem is due to Ismail and comprises in an abstract form the fundamental 
fact behind the theory of basic analogs of Fourier series on a g-quadratic grid [5] . 

Theorem A 15] The system {F n (x)} is orthogonal and complete in L 2 (li). 

To give an idea of what is involved in this statement, we sketch Ismail's argu- 
ment. Since, by (|1.3p . rk(t n ) are the Fourier coefficients of F n in the basis {ukPk}, 
the use of Parseval's formula gives 



/oo 
F n (x)F m (x)dLi(x) = y^ y r k (t n )r k (t m ) 



P(Xn) 

and the orthogonality relation is proved. To show the completeness, choose / € 
L 2 (li) and assume J F n (x) j '(x)d/j,(x) — Oforalln. Again Parseval's formula implies 
SfcLo fkUk^kitm) = for all m, where are the Fourier coefficients of / in the 
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basis {pk}- Now the completeness of {ru} implies fk = 0. Therefore / = almost 
everywhere in L 2 (/i) . 

In [16] (see also section 24.2 of the monograph |21j ), Ismail posed the problem 
of studying the mapping properties of operators with kernels defined as above and 
conjectured that there was a reproducing kernel Hilbert space structure behind 
these operators. We will show that Ismail's conjecture is true. Our approach will 
reveal a reproducing kernel structure reminiscent of the well known structure of the 
Paley- Wiener space of functions bandlimited to a real interval. However, even in 
the case when the system {F n (x)} is the set of the complex exponentials, we obtain 
results that, as far as our knowledge goes, seem to be new. When the system 
{F n (x)} is the set of basis functions of the g-Fourier series constructed with the 
function £ q (x;it), we will obtain results that complement the investigations done 
in [29] and [19] . In particular it will be shown that a sampling theorem related 
to the one derived in [19] lives in a reproducing kernel Hilbert space and that the 
correspondent g-analogues of the Sine function provide an orthogonal basis for that 
space. 

Before outlining the paper we want to make clear that the techniques we are 
using already exist in some antecedent form. In particular, section two contains 
ideas from Ismail theory of generalized q- Fourier series and from Saitoh 's theory of 
linear transformations in Hilbert space. However, their particular combination here 
leads to new conclusions and sheds new light in the emerging theory of g-Fourier 
series. It reveals an elementary structure underlying many systems involving several 
special functions simultaneously. The results in section 2.2 and section 3, 4 and 5 
were never explicitly stated before. 

The paper can be summarized as follows. The next section contains the de- 
scription of the reproducing kernel structure behind the abstract setting of theorem 
A. An integral transformation between two Hilbert spaces is defined in the context 
of Saitoh 's theory of linear transformations, basis for both spaces are provided and 
the formula for the reproducing kernel of the image Hilbert space is deduced. In 
this context a sampling theorem appears, generalizing the one in [19] . The remain- 
ing sections consider three applications of these results, using specific systems of 
orthogonal polynomials as well as Bessel functions and their generalizations. The 
first application is associated to formula and systems of complex exponentials. 
The second application is linked to (II. 2ft and to systems of curly g-exponentials, 
and we write the reproducing kernel as a 2<t>2 basic hypergeometric function. These 
two examples explore the interplay between Lommel polynomials and Bessel func- 
tions and the corresponding relations between their g-analogues. In the last section 
we consider a construction of a general character, designed originally in the papers 
[20] . [17] and [15] . It allows to extend the interplay between Bessel functions and 
Lommel polynomials to a more general class of functions. Using this construc- 
tion we will make a brief discussion about the application of our results to spaces 
weighted by Jacobi weights and their g-analogues. 

2. The reproducing kernel structure 

Let if be a class of complex valued functions, defined in a set X C C, such 
that if is a Hilbert space with the norm of L 2 (X,fi). The function k(s,x) is a 
reproducing kernel of H if 

i) k (., x) £ H for every x G X; 
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ii) f (x) = (f (.) , k (., x)) for every f £H,xeX. 

The space H is said to be a reproducing kernel Hilbert space if it contains a 
reproducing kernel. From a structural point of view, the correct approach to the 
study of our problem is via Saitoh's theory of linear transforms of Hilbert space. 

2.1. Preliminaries on Saitoh's theory of linear transformations in 
Hilbert space. This theory can be found in works by Saitoh [26], [27j and we 

proceed to give a brief account of the results that we are going to use. An account 
of the results quoted in this section can also be found in Higgins recent survey [13j . 

For each t belonging to a domain D, let K t belong to H (a separable Hilbert 
space). Then, 

k(t,a) = (K t ,K.) H 

is defined on D x D and is called the kernel function of the map K t . Now consider 
the set of images of H by the transformation 

(Kg)(t) = (g,K t ) H = f(t) 

and denote this set of images by R{K). The following theorem can be found in 
[21]: 

Theorem B The kernel k(t, s) determines uniquely a reproducing kernel Hilbert 
space for which it is the reproducing kernel. This reproducing kernel Hilbert space 
is precisely R(K) and it can have no other reproducing kernel. 

Now, suppose that {K t }(t £ D) is complete in H so that K is one to one. Then 
we have 

\\ K 9\\r{ K ) = h\\ H - 

The following theorem, due to Higgins [12] . will be critical on the remainder. We 
will use only the following " orthogonal basis case" , a special case of the theorem in 

m- 

Theorem C With the notations established earlier, we have: If there exists 
{t n } (n £ I C Z) such that {K tn } is an orthogonal basis, we then have the sampling 
expansion 

in R(K), pointwise over I, and uniformly over any compact subset of D for which 
\\K t \\ is bounded. 

2.2. The reproducing kernel for (/-Fourier type systems. In this section 
we will show the existence of a reproducing kernel structure behind the abstract 
setting of Theorem A. The results will follow from the study of the mapping proper- 
ties of an integral transform whose kernel is obtained from the sequence of functions 
{r k } and {p k }. 

Theorem 1. There exists a kernel K(x, t) satisfying the requirements of Saitoh 's 
theory of linear transformations, such that K(x,x n ) = X n F n (x), where {F n (x)} is 
the orthogonal sequence of functions in Theorem A and {X n } is a sequence of real 
numbers. 

Proof. Our first technical problem comes from the fact that, when {r k } is 
a discrete system of orthogonal polynomials with a determinate moment problem, 
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then {rk(t)} £ I 2 if and only if x is a mass point for the measure of orthogonality. 
For this reason the series 



fc=0 



rk (t)p k (x)u k 



would diverge if t is not such a point (this is pointed out in Section 5 of |15j ). 
Since we want our kernel to be defined for every t, we will assume the existence of 
an auxiliary system of functions {Jk{x)} € I 2 for every t real, and such that every 
function J k interpolates r k at the mass points {x n } in the sense that 

(2.1) J k (x n ) = A„r fe (— ) 

for every k = 0, 1, ... and n = 0, 1, ...and some constant A n independent of k. Now 
we can use the functions Jk{t) to define a kernel K(x,t) as 

oo 

(2.2) X(M) = J! Jk{t)Pk(x)u k . 

k=Q 

Such a kernel is well defined and belongs to L 2 (p), since it is a sum of basis functions 
of L 2 (ll). From dOJ), flSj]) and we have 

OO 

i*r(x,£ n ) = ^2 Jk{x n )pk{x)u k 
k=0 



oo ^ 

A„^r fc ( — )p fc (a 



k=0 ™ 
= A„i^(x). 

□ 

Remark 1. Observe that theorem 1 and theorem A with t„ = — show that 
K(x,x n ) is an orthogonal basis for the space L 2 (fi). 

Remark 2. In the abstract formulation it may not be clear why the constant X n 
must be present. Actually the construction would work without it, but for technical 
reasons that will become evident upon consideration of examples we prefer to use it. 
Otherwise, careful bookkeeping of the normalization constants would be required in 
the remaining sections. 

Remark 3. It will be seen in the last section that a general constructive method 
is available in order to find the function J k under very natural requirements on the 
polynomials r k - 

Now define an integral transformation F by setting 

(Ff)(t) = J f{x)K(x,t)dn(x). 

We will study this transform as a map whose domain is the Hubert space L 2 (p). 
Endowing the range of F with the inner product 

(2-3) (Ff,Fg) F{L2M) = (f,g) L2M 

then F(L 2 (/j,)) becomes a Hilbert space isometrically isomorphic to L 2 (fi) under 
the isomorphism F. Using Saitoh's theory with D — R, Kt — K(., t), H — L 2 (/j,) 
and R{K) = F(L 2 (^i)), we obtain at once: 
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Theorem 2. The transform F is a Hilbert space isomorphism mapping the 
space L 2 (/j,) into F(L 2 (fi)). The space F(L 2 {p)) is a Hilbert space with reproducing 
kernel given by 



(2.4) k(t,s) = J K(x,t)K(x,s)dfj,(x). 

An interesting feature of this particular setting is that every function in F(L 2 (p,)) 
has two different expansions: One is the sampling expansion naturally associated 
with the reproducing kernel structure, the other is the expansion in the basis J n (x). 
It should be remarked that, in most of the previously known sampling expansions, 
these two expansions were the same. These expansion results are summarized in 
the next theorem. 

Theorem 3. Every junction f of the form 

(2.5) f{x) = J u(t)K(t,x)dti(t) 
with u £ L (fi), admits an expansion 

(2.6) /(*) = > - 
where the coefficients are given by 



ri=0 



U n (u,p„(.)) L 2 (p) 



and a sampling expansion 

(2-7) /(a)= £ /(tn )^> 

The sum in \2. 7[ ) converges absolutely. Furthermore, it converges uniformly in 
every set such that H-KX-j^ll/^f^) * s finite. 

Proof. We already know by default that {p n (x)} is a basis for L 2 {^l). It 
remains to prove that {J n {t)} is a basis for F(L 2 (p)). Observe that 



(Fp n )(t) = Pn (x)K(x,t)dn(x) 



k=0 



Jkit)u k I p n (x)p k (x)d[l(x) 



— ^7n(^)^n- 

Since {p n {x)} is a basis for L 2 (/i) and F is an isomorphism between L 2 (p) and 
F(L 2 (/i)), then {u n J n (x)} is a basis for F(L 2 (/i)). To prove the last assertion of 
the theorem, observe that the function / defined by (|2.5|) belongs to F(L 2 (/j,)) and 
therefore can be expanded in the basis {u n J^ n (x)} . The Fourier coefficients of this 
expansion are 

a n = (f,u n J n (.)) F{L 2 M) = (Fu,F(p n (.))) F{L2{tl)) = (u,p n (.)) L 2 M 

where we have used (|2.3|) in the last identity. The sampling expansion follows from 
applying theorem C to our setting and using remark 1. □ 
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Remark 4. The construction of this section has never appeared before in the 
literature, but it is reminiscent of the reproducing kernel structure of the Paley- 
Wiener space. In the classical situations (see for example [8] and [24] for an account 
of these constructions with several examples ) generalizing this structure, there is an 
integral transform whose kernel is defined as 

(2.8) K{x,t) =^S fe (t)e fe (x) 

where ek{x) is an orthogonal basis for the domain Hilbert space and Sk(t) is a 
seguence of functions such that there exists a sequence {t n } satisfying the sampling 
property 

(2.9) Sk(t n ) = a n S n ^ k 

As an instance, take Sk(t) — an d &k{x) — e lkx . Then 



e itx _ \^ sin n(t — k) ^ ikx 



£c «■(*-*) 



and K(x, t) is the kernel of the Fourier transform. The corresponding reproduc- 
ing kernel Hilbert space is the Paley- Wiener space. The root of these ideas is in 
Hardy 's groundbreaking paper [11] . For an application of this classical set up to 
Jackson q-integral transforms and the third Jackson q-Bessel function, see pQ. In 
our construction we made a modification of this classical setting: Instead of the se- 
quence of functions Sk, with the sampling property i2.9\) . we considered a sequence 
of functions {j7fc}, interpolating an orthogonal system {r^} in the sense of \2.1\) . 
And we have seen that the essential properties of classical reproducing kernel set- 
tings are kept. However, this modification allows to recognize a class of reproducing 
kernel Hilbert spaces that were obscured until now. This will become clear in the 
next section. 

3. The Fourier system with ultraspherical weights 

The nth ultraspherical (or Gegenbauer) polynomial of order v is denoted by 
C^{x). These polynomials satisfy the orthogonality relation 

L n {x)L m {x){l x) dx- nl{v + n)v{y) <W 

and form a complete sequence in the Hilbert space L 2 [{— 1, 1), (1 — a; 2 )" -1 ' 2 ]. For 
typographical convenience we will introduce the following notation for this Hilbert 
space: 

H v = L 2 [{-1,1),(1 - x 2 )"- 1 ' 2 ]. 
The Bcssel function of order v, J v {x), is defined by the power series expansion 



JS, (-Y) n / z\ v + 2n 

<«> M*) = Y. n!r l + l +1} (|) ■ 

n— 

The nth Lommel polynomial of order v, denoted by h njU (x), is related to the Bessel 
functions by the relation 

(3.2) Ju+k(x) = hk,u{-)Ju(x) - hk-i, v -i(-)Ju-i(x). 

x x 
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The Lommel polynomials satisfy the discrete orthogonality relation 

00 1 1 1 1 

S ^^'^l^^'^^ = 2{u + n + l) 6nm 
and the dual orthogonality 

00 11 

2J 2(l/ + n + l)h klV+ i(±- )hk,u+l(±- ) = (jv,k) 2 8nm- 

3v,n jL/ : m 

They form a complete orthogonal system in the l 2 space weighted by the discrete 
measure with respect to which they are orthogonal. We will use these two complete 
orthogonal systems in our first illustration of the general results. Set 



Pfc(x) = V^)^ k[X) 



and 



r-fc(t) = \j2(v + n)hk,v-i(t). 

Consider also 

J k (t) = y/2(v + n)J v+k (t). 

Denote by j v ± the kth zero of the Bessel function of order v. Substituting x = j u . 
the following interpolating property is obtained 



(A^J V-Hj ) - — j p \- 

The interpolating property (|3.3[) will play the role of (|2.1j) with A„ = — j — 4^ — r. 
Consider also the sequence of complex numbers {u n } defined as 

Uk = i k 

and set 



Now, theorem 1 tells us that K"(x,j Uin ) is an orthogonal basis of the space H v ' . 
Moreover, formula (jl.ip implies 

f 

and therefore we can think of T(y) (|) K v (x, t) as a one parameter generalization 
of the complex exponential kernel that may be worth of further study. In the 
special case v = h we recover the well known orthogonality and completeness of 
the complex exponentials {e l7rnx } in L 2 (— 1, 1). 

The transformation F is defined, for every / E H u , as 

(Ff)(t) = J f(x)K»(x, t)(l - i 2 )" -172 ** 

and theorem 2 gives that the reproducing kernel of F(H V ) is 

k"(t,s)= [ K v {x,t)K»{x,s){l - x 2 ) v - l ^dx. 
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When v = 1/2 this becomes 

iryts sin(i — s) 



k*(t, s) 



(t-s) 



Since {y ^^y^-C^(t)} forms a basis of the space H u , then also {i n yj2{v + n) J v+n {x)} = 
F {\J ^j0y~Cn{t)} is a basis of the space F (H lJ ). In this situation theorem 3 reads: 
Theorem 4. Let f be a function of the form 

(3.4) f(t) = ^ u(x)K u (x, t)(l - x 2 f- 1/2 dx 
where u G H v . Then f can be written as 

oo 

(3.5) f(t) = J2 a nJ»+n(t) 

n=0 

with the coefficients a n given by 



(3.6) a n = i n {u + n)J^- f u{x)C»{x){l x^^dx 



i 



and also as the sampling formula 

k u {t, ju,n) 



/(*) = £ f^u,n): 



Remark 5. Expansions of the type VS. 5)) are known as Neumann series of 
Bessel functions (see chapter 16 of 25 ). 

Remark 6. When v = 1/2 the above sampling theorem states that every func- 
tion of the form 

b > I ' " ixt. 



fit) = ( 2 ) / u( x y xt dx, 



with u G L 2 [{— 1, 1)], can 6e represented as 

„. . „,„ „ / t sm(t — 2im) 

/(*) = g/(2^— { l 2m) K 

4. The a-Fourier system with g-ultraspherical weights 

We proceed to describe the g-analogue of the previous situation. Choose a 
number q such that < q < 1. The now classical notational conventions from [9] 
and |21j for o-infinitc products and basic hypergeometric series will be used often. 

The q-exponential function that we talked about in the introduction is defined 
in terms of basic hypergeometric series as 

(-^;g^ )oo x ( qie i0 ,qie~ ie 
iqt 2 ;q 2 



£qfat) = —7—7— 2'.' i 



q*,-t 



where x = cos 0. The continuous q-ultraspherical polynomials of order v are denoted 
by C%(x; q"\q) and satisfy the orthogonality 

l^(x;q v \ q )C^(x ]q v \ q )wix; q ^ \ q)dx = ^Z^f^J ^n 
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where the weight function w(x; (3 \ q) is 

uj(cos#; /? o) = - — — — - — 7- „. a - — ^— — - — , (0 < < 7r). 

The polynomials {C^(x;q v \q)} form a basis of the Hilbert space defined as 

HZ = L 2 [(-i,i),w(x;q v \q)}- 

The second Jackson g-Bessel function of order v is defined by the power series 

^ [X ' q >- (q-^U 2J l) (q;qUq^;q) n q 

Since this is the only g-Bessel function to be used in the text, we will drop the 

(2) 

superscript for shortness of the notation and write J v (x\q) = (x;q). The q- 
Lommel polynomials associated to the Jackson g-Bessel function of order v are 
denoted by h niV —x(x; q). These polynomials were defined in |14j by means of the 
relation 

(4.1) q ^+n{n-X)/2 j ( q) = h ( 1 ^ j v ^.. j _ K _ (I. q ) J v _^ x . g ). 

X X 

The q-Lommcl polynomials satisfy the orthogonality relation, for a certain constant 
A n (v + 1) which is not explicitly known (see section 14.4 of |21j): 

°° A n (v+1) 1 1 qav+n(n+l)/2 

/ , — TT^h n v+ i{±- —■ 1 q)h m , v+ i(±- 7-r',q) = — , , -, S 



win 



and the dual orthogonality 

- (l_^+i) 1 1 _ (i,,„(g)) 2 

2-, 9 », + n(n + l)/2 ^,,+U± ^ , <Zjft m ,„ + lt±^ , <Z) - ^ + 1} 0„ m . 

Consider 



{l-g k+v )(q;q)k „ 

ft(l) lR» Cfc(li9k) ' 



and 



Jfc(t) = 

The parameters will be given by 



Denote by j v ,k{q) the fcift. zero of J„(x; g). Setting t = j»,k(q) in ()4.1j) we have the 
interpolating property 

(4.2) hk i/ _ l( -J—. q) = -q*>>+Hk-i)/2 J^Aq)', <?) 



REPRODUCING KERNELS 



11 



This means that in Ij2.1j) we must take A„ = — -j — ( q y q y ln this context, the 
kernel K(x, t) is given as 

+\ \^ jkMiv+k-p/i 0- ~ q k+v ) / (g; g)fc T /„,, v r / v. \ 
* g (aj,i) -2^« 9 (1 - g") V (g 2t/ ; g) fc +fc( g)Cfc(x ' g l<z) 

and the use of (|1.2|) gives, when v =\, 

(4.3) tf| (x, t) = (-g^ 2 ^ g^)oo(g" +1 ; g)oo ^ g , it) _ 

(9; 9)00 

The basis functions of the domain space are 

F n (x) = K^x,j u , k (q)). 

When v = \ this gives that {£ q (x;iji n _ 1 (q))} is orthogonal and complete in H%. 
This is the case of the g-Fourier series studied in [29] . Now define the transform 

(4.4) (F;/)(<) = J f(x)K»(x,t)w(x;q" \ q)dx 
for every / £ H q . We have 



l V mt 9 4 ■ Wfc(2 * 5?)=F - Vvf 1 -^")* Cfe( * ;(H<?) J 

and { J iy+n (t; q)}^ =0 is a basis of the space (-ff,). We can also state the q- 
analoguc of the expansion in theorem 5 and obtain a q- Neumann expansion theorem 
in (j-Bessel functions. 

Theorem 5. Let f be a function of the form 
(4.5) /(*)=/ u{x)K u Jx,t)w{x;q u \ q)dx 



where u G . Then f can be written as 

q j 

00 

/(*) = ^2a n J v+n (t;q) 

n=0 

with the coefficients a n given by 



(l-q k+v ) / (q;q) k f 1 „ 

' u(x)C n (x;q \q)w{x;q \q)dx. 



(i-g") V (Q 2v ;<!hJ-i 

Once more we know by theorem 2 that the space F q [H q j is a space with 
reproducing kernel k q (t, s), given by 

= y K^(x,tjK g 7 {^s)w(x;q u \ q)dx. 
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When v — i we can compute the reproducing kernel in an explicit form using the 
following integral, evaluated in [18j : 



(4.6) / £ 9 (cos0;a)£ g (cos0;/3) 



( 7e 2ie ;7e -2i8. g ) c 



d9 



2v('y,qj,-af3q 1/2 ;q) 00 ^ 
(q,l 2 ;q)oo(qa 2 ,qf3 2 ;q 2 )oo 



^ 2 a/(3,-q^ 2 P/a 
qj, , -a^g 1 / 2 



q, -af3-yq 1/2 



Theorem 6. The space F q 2 (^Hq^j * s a space with reproducing kernel kq (t,s), 
given by 

V /2 ;?) 3 oo /2 "" 



(4.7) 



x {q 1 l 2 ,-tsq 1 l 2 ; q) x (ts)h<t> 2 



(q;q)l 



q y 2 t/s,~ q y 2 s/t 

7 3/2 



q, -tsq 



q"'~,-tsq 

PROOF. Applying theorem 2 and (|4.3[) we know that kg (t, s) is given by 



(q v+1 ;q)c 



l 2 



(-gt ,-qs ;q ) 00 (ts) 



x y £q(x; it)£ q (x; — is)w(x; q v \ q)dx. 

Make the substitutions x — cos0, it — a, is — f3, and q v = 7 in (|4.6I) . Then (|4.7p 
follows. □ 

As in the preceding sections we can formulate a sampling theorem. 

Theorem 7. Every function of the form |^.5| ) admits the expansion 



(4.8) 

w/iere t fc = i^i. 



Remark 7. XTie functions k i^ t ' k tk ^ are related to the functions Sinc q (t,k) 

from |19) . In [19] is shown that there is a sampling theorem associated to the 
Sinc q (t, k) functions that can be written as an interpolating formula of the Lagrange 
type. The above discussion shows what is the sampling theorem associated to the 
reproducing kernel Hilbert space and its relation to the one obtained in 



Remark 8. By the proof of Theorem B (see (13] ), the functions fc^(t fc 't fc ) are 

orthogonal in the space F q 2 (Hq ). This result is a q-analogue of the important fact, 
proved by Hardy in [llj . that the functions ~7^~fc^ are orthogonal in the classical 
Paley-Wiener space. 

Remark 9. Important information concerning the zeros of the second Jackson 
q-Bessel function, that appear as sampling nodes in the expansion 14- 8\) , was ob- 
tained very recently by Walter Hayman in [10] using a method due to Bergweiler 
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and Haymann [4] : the asymptotic expansion 

n 

f v , k (q) = 4q 1 - 2n -"{l + Y, b ^ +° 



Jn+l)k 



} 



as k — > oo, with the constants h v depending on a and q, holds. Therefore, for big 
k, the sampling nodes are exponentially separated in a similar way to what was 
observed in pQ and [2]. In the case where v — ±j, the zeros were studied by Suslov 



5. A generalization 

We begin this last section describing a formal approach generalizing the situ- 
ations studied in the two previous sections. This formal approach was initiated in 
[20j and [17] with the purpose of finding functions to play the role of the Lommcl 
polynomials in more general situations, and was studied further in [15] . In the 
context studied in this paper, it will be of particular relevance, since it gives a con- 
structive method to find the functions Jk satisfying (|2.1[) . Let {f n ,v\ be a sequence 
of polynomials defined recursively by fo,v{x) — 1, fx,v{x) = xB v and 

fn+l,u(x) = [xB n+u ]f n>v (x) - C n+u -if n -l,v(x). 

Assuming the positivity condition B n+l ,B n+v+ iC n + v > (n > 0) and the conver- 
gence of the series X)^o b + b + " + ~ ' ^ can ^ e snown j using facts from the general 
theory of orthogonal polynomials, that the polynomials /„ „ are orthogonal with 
respect to a compact supported discrete measure and that the support points of 
this measure are , where the x n v are the zeros of an entire function J satisfying 

(5.1) C v ...C v +n-iJ(x; v + n) — J(x\ v)fn,v{~) - J(x; v - l)/„_i v+1 (-). 

x x 

The dual orthogonality relation of the polynomials f n ,v{x) is 

5 > o\ 7 i 7\fn,v+x( )fn,v+l( ) = ~7~, ""TT^j.fc 

^2A„(z/ + l) x v .k x u j Ajiy + l) 

for some constants Aj{v + 1) and \ n (v + 1). (for the evaluation of these constants, 
as well as other parts of the argument missed in this brief sketch, we recommend 
the reading of section 4 of [15] ). From (]5.ip and the above analysis we can obtain 
the interpolation property 

J \X n ,u\ V + k) = — Jk-l, v +\[ 



Cv"-Ci/-{-n—l x nu 
Therefore, in the language of the second section we can set 

-B k 



Jkit) = y/\^u)-^J(t, k) 

r k (t) = y/J^)^Lf niU (t) 
. _ -J(x n y, v - f) 

n n n 

^v-^u+n— 1 
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and define the kernel 



K(x,t) 



^ UkV^k(v) ^ J(t, k)p k (x) 




where \uk\ = 1 and {p n (x)} is an arbitrary complete orthonormal system in L 2 (n). 
As before, the kernel K (x, t) can be used to define an integral transformation be- 
tween two Hilbert spaces. We could now apply the machinery of section 2 and 
provide a reproducing kernel structure and a sampling theorem by means of an 
integral transform with the above kernel. However, no simplification would occur 
on the absence of proper addition formulas for the kernel K (x, t). Choosing families 
of orthogonal polynomials f n ,v{t) and p n {x) in a way that such addition formulas 
exist is the topic of the second problem in [16j . 

Operators weighted by the Jacobi weight function (1 — x) a (l + x) 13 on the 
interval [—1, 1] can be studied by using the Jacobi polynomials and the Wimp 
polynomials in a similar fashion to what was done in section 4 of |20[ section 4]. 
Reasoning as before, we can use the resulting formulas to generalize the results 
in the third section to Fourier systems with Jacobi weights, although no major 
simplification seems to occur. Results very similar to those of section 3 would 
follow, with an extra parameter; expansions in series of \F\ replace the Neumann 
expansions and sampling theorems with sampling points located at the zeros of 
these ii<i can also be derived. 

In section 6 of [IT] . formula (6.13) is a generalization of (|1.2[) . This formula 
involves continuous q-analogues of the Jacobi polynomials defined via the Askey 
Wilson polynomials 3 and a g-exponential function with an extra variable. This 
more general and complicated situation should be studied in a future work once 
more summation formulas are known. Also, following Ismail's suggestion in Prob- 
lem 24.2.1 of |21j may lead to the discovery of other examples of structures fitting 
to the setting described in this paper. 
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